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SUMMARY 


A  very  almple  algorithm  for  finaing  a  maximal  flow 
and  minimal  cut  in  a  transportation  nctwora  la  described 
it  la  then  applied  to  obtain  an  efficient  computational 
routine  for  the  Hitchcock  distribution  problem. 
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A  SIMPLE  ALGORITHM  FOR  FINDING  MAXIMAL  NETWORK  FLOWS  AND 
AN  APPLICATION  TO  THE  HITCHCOCK  PROBLEM 


L .  R  .  Ford  ,  Jr  . 
D.  R.  Fulkerson 


INTRODUCTION 

The  network— flow  problem,  originally  posed  by  T.  Harris 
of  The  RAND  Corporation,  has  been  discussed  from  various  view¬ 
points  in  [l]  ,  [d]  ,  [  }1  ,  '  ,  ,  ,  .  .  -  >  arises  naturally 

in  the  study  cf  transports  1  i  on  natvoiva;  it  may  be  stated  in 
the  following  way.  One  is  given  a  otwork  of  directed  circa 
and  nodes  with  two  distinguished  nados,  called  source  and  3 Ink , 
respectively.1  All  other  nodes  are  called  intermediate .  Each 
directed  arc  in  the  network  has  associated  with  it  a  rionnegatlve 
integer,  its  flow  capacity.  Source  arcs  may  bo  assumed  to  be 
d lrectod  away  from  the  source,  sink  arcs  into  the  sink.  Subject 
to  the  conditions  that  the  flow  in  an  arc  is  in  the  direction 
of  the  arc  and  uoes  not  exceed  its  capacity,  and  that  the  total 
flow  into  any  Intermediate  node  is  equal  to  the  flow  out  of  it. 

It  la  desired  to  find  a  maximal  flow  from  source  to  3ink  in  the 
network,  l.e.,  a  flow  which  maximizes  the  sum  of  the  flows  in 
source  (or  sink)  arcs. 

1  A  problem  in  which  there  are  several  sources  and  sinks, 
with  flows  permitted  from  any  source  to  any  sink,  Is  reducible 
to  a  single — source,  single— sink  problem. 


For  example ,  consider  the  networx  ol  Pig  •  1  witn 
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30urce  P1#  oiriK  P4,  ami  arc  capacltic3  aa  indicated 

let  x,  <  denote  tne  flow  from  P,  to  P  ,  ,  the  probi’M.  13 
J 

imize  x  1 2  4-  x  1 3 ,  the  total  flow  leaving  Pw  sub,  ect  t< 
equationa  ajid  inequalities 

X)2  4  Xaa  “  X  2  3  *2*  , 

Xt)  +  Xgj  *  X32  4-  Xjt  , 

O  <  xl2  <  1, 

0  h  x  j  3  <  1  , 

0  <,  x  2  3  <  1 , 

0  £  x»2  P  2  > 

0  £  *24  <  1» 

0  <;  x34  <  >. 

In  general,  if  we  let  Pi  be  the  source,  P^  tne  sln^, 
required  to  find  x1J(l,J  -  !.-••«)  which  maximize 


If  we 
to  max— 
the 


we  are 


CO 


subject  to 

(2)  ^  ““  *jj_)  "  ®  (^  "  2,...,n—  1)  , 

0  ^  *ij  -  °1J  * 

where  the  are  given  nonnogative  integers,  and.  In  particular, 

cil  ~  cnJ  *  0  for  a11  ^  * 

TYiia  la  of  course  a  linear  programming  problem,  and  hence 
may  be  solved  by  Dantzig's  simplex  algorithm.  In  fact,  the 
simplex  computation  for  a  problem  of  this  kind  ia  particularly 
efficient,  since  it  can  be  shown  that  the  sets  of  equations  one 
■olvea  in  the  process  are  always  triangular  [i] .  However,  for 
the  flow  problem,  we  shall  describe  what  appears  to  be  a  con¬ 
siderably  more  efficient  algorithm;  it  is,  moreover,  readily 
learned  by  a  person  with  no  special  training,  and  may  easily 
be  mechanized  for  handling  large  networks .  We  believe  that 
problems  involving  more  than  t?00  node3  and  **,000  arcs  are 
within  reach  of  present  computing  machines  . 

Of  some  theoretical  interest  la  the  fact  that  the  procedure 
assures  one  of  obtaining  a  strict  increase  in  the  total  flow  at 
each  step  (in  contrast  with  the  simplex  method).  In  addition, 
the  Hltchoook  problem  can  be  solved  via.  the  flow  algorithm  in  a 
way  whioh  naturally  generalizes  the  combinatorial  method  recently 

t 

proposed  by  Kuhn  [7]  for  the  optimal- assignment  problem.  Informal 
tests  by  hand  indicate  that  this  way  of  solving  the  Hitchcock 
problem  is  extremely  efficient. 
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*  *  COMPUTATION  OF  MAXIMAL  PLOW,  AN  EXAMPT  R 

In  order  to  illustrate  the  computational  procedure,  before 
describing  it  in  general  terms,  let  us  return  to  the  example  of 
Fig.  1.  Start  out  with  any  flow  from  source  to  sink,  say  a  flow 

of  1  along  the  chain  Pt  P*  pa  p4.  ye  describe  this  flow  by  Fig. 
2  below. 


© 


Fig .  2 

where  the  remaining  or  unused  capacity  from  to  PJ  ia  denoted 
by  the  number  slosest  to  P1  on  the  arc  P±  P^  .  Notice,  for  example, 
that  the  new  capacity  from  P.  to  Pa  has  been  decreased  by  1 ,  and 
the  capacity  from  P,  to  P,  has  been  increased  by  1.  That  is  to 
say,  we  must  new  allow  the  possibility,  in  constructing  successive 
flows,  of  imposing  a  flow  of  2  from  P8  to  P,,  one  unit  of  which 
would  cancel  the  present  now  of  1  from  ps  to  Pa.  Now  start 
with  the  source  and  look  for  nodes  which  may  be  reached  "in  one 
step  by  arcs  of  strictly  positive  remaining  capacity.  Here  we 
may  proceed  only  to  P*.  When  a  node  has  been  reached,  label  it 
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in  some  fashion;  then  take  a  labeled  node  not  previously 
examined,  and  look  for  unlabeled  nodes  which  may  be  reached 
in  one  step  by  arcs  of  positive  remaining  capacity.  -Here  we 
may  proceed  from  P*  only  to  P».  Repeat  the  procedure  until 
either  (a)  the  sink  has  been  labeled,  or  (b)  no  further  nodes 
may  be  labeled  and  the  sink  has  not  been  labeled.  In  our  example, 
case  (a)  occurs  at  the  next  step.  We  may  now  search  back*  through 
labeled  nodes,  locating  a  chain  from  source  to  sink  along  which 
an  additional  flow  may  be  imposed.  Hare  we  have  the  chain  PxPaPaP 
and  a  flow  of  2  may  be  imposed.  We  thus  obtain  Pig.  5  b«low 


Pig.  2 


from  71g.  2  by  alternately  subtracting  and  adding  2  to  the 
numbers  encountered  along  the  chain  PiPb^sP*.  Next  repeat  the 
labeling  procedure.  This  time  we  can  label  only  Pi  and  Pa, 
and  are  thus  in  case  (b).  T*ia  means  that  we  have  obtained  a 
maximal  flow,  depicted  in  Pig.  • 


*  In  the  general  description  of  the  algorithm  to  be  given 
presently,  we  will  carry  enough  information  along  to  make  the 
backward  search  unnecessary . 


>3 
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Pig.  4 

The  proof  that  this  flow  is  maximal  in  this  case  is  • 
immediate ,  for  observe  that  the  directed  arcB  PtP»,  PaPs,  Pa?4 
leading  from  labeled  nodes  to  unlabeled  nodea  aro  all  saturated 
and  form  a  cut  in  the  network;  l.e.,  every  directed  chain  from 
source  to  sink  contalna  one  of  these  dlreoted  arcs.*  Slnoe  it 
la  clear  that  the  sum  of  capacities  of  arcs  forming  a  cut,  which 
we  refer  to  as  the  value  of  the  cut,  is  an  upper  bound  for  flow 
values,  and  since  we  have  achieved  equality  of  flow  value  and  a 
cut  value,  the  flow  is  maximal  and  the  directed  arcs  Pi  Pa,  PaP», 
PaP4  constitute  a  minimal  out,  i.e.,  a  cut  of  minimal  value. 

2 .  THE  MINIMAL  CUT  THEOREM 

A  none one truotive  proof  of  the  minimal  cut  theorem,  which 
asserts  equality  of  maximal  flow  value  and  minimal  cut  value, 
has  been  given  by  the  present  writers  in  [3] .  Subsequently  a 
constructive  proof  baBed  on  the  simplex  criterion  of  linear 

*Tho  definition  of  cut  given  here  corresponds  to  that  of 
disconnecting  Bet  given  in  [5] • 
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programming  was  developed  [2,5>  The  algorithm  which  we  describe 
more  formally  In  the  next  section  aiao  provides  a  constructive 
proof  of  the  theorem.  Like  the  simplex  algorithm.  It  produces 
not  only  a  maximal  flow  but  a  minimal  cut  as  well.  This  will 
be  Important  for  our  application  to  the  Hitchcock  problem. 

It  should  perhaps  be  pointed  out  that  an  undirected  problem, 
by  which  we  mean  that  the  directions  of  flow  in  intermediate 
arcs  are  not  specified,  so  that  the  capacity  constraints  on 


these  arcs  are  of  the  form 


O)  *u  *  *ji  *  cij  J 

presents  nothing  new,  Bince  we  may  replace 
by  a  pair  of  oppositely  directed  arcs,  each 
to  that  of  the  original  arc;  i.e.,  replace 


• , n — 1 ,  1<J )  , 

aach  undirooted  arc 
with  capacity  equal 
(3)  by 


X1J  $  C1J  ' 

XJ1  *  C1J  ’ 

For,  given  X*  -  (x'j)  satisfying  (4)  and  the  conservation  equations 
at  intermediate  nodes,  setting 

(5)  -  max  <xij  -  xji,0) 

yields  an  equivalent  flow  X  -  (x^)  aatl.fylng  theee  equatlona 
and  (3)  . 

The  minimal  cut  theorem  is  true  for  undirected  networks  aa 
well  as  for  directed  networks,  or,  more  generally,  for  mixed- 
networks  in  which  some  Intermediate  arcs  are  directed,  others  not, 
the  obvious  changes  in  definitions  of  cuts  and  chains  having  been 


made.  This  follows  from  the  comments  above  and  the  fact, 
easily  proved,  that  the  minimal  cut  value  is  the  same  for  a 
mixed  network  and  its  equivalent  directed  network.  The  theorem 
is  still  valid  when  oapacity  constraints  on  nodes  are  admitted, 
where  a  cut  now,  of  course,  Includes  nodeB  as  well  as  aros. 

This  may  be  proved  by  splitting  each  node  into  two  nodes  as 
suggested  In  Fig.  5 


Fig.  5 


and  placing  the  capacity  o  of  the  old  node  on  the  new  directed 
arc  Joining  the  two  new  nodes,  thus  obtaining  an  equivalent 
network  with  capacity  constraints  on  arcs  only.  For  a  direct 
proof  of  the  theorem  in  this  general  form,  see  [2]  . 

TKB  ALGORITHM 

We  start  the  computation  from  any  convenient  initial  flow 
whatsoever.  The  initial  flow  is  used  to  define  a  starting 
matrix  A  -  (a^)  by  letting  be  the  capacity  of  the  arc  from 
Pi  to  Pj  diminished  by  the  flow  from  P^  to  P^  and  increased  by 
the  flow  from  Pj  to  P^.  If  no  other  flow  is  readily  available, 
one  may  st^rt  with  the  zero  flow,  corresponding  to  A  ■  C,  where 
c^j  is  the  original  oapacity  from  P^  to  . 
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(It  la  the  opinion  of  the  authora  that  If  the  problem 
la  given  in  matrix  form,  no  special  attempts  should  be  made 
to  obtain  a  good  starting  solution.  If,  on  the  other  hand, 
the  problem  can  be  plotured  readily  as  a  linear  graph,  a  "flooding" 
idea  described  in  [l]  might  be  used  to  obtain  a  starting  flow. 

By  following  the  approach  suggested  in  [l] ,  which,  however, 
calls  for  the  exercise  of  Judgment,  an  initial  flow  can  be 
obtained  that  often  is  optimal  for  simple  networks.  If  not, 
it  might  be  used  as  a  good  starting  point  for  initiating  the 
procedure  given  in  this  paper.) 

We  are  assuming  that  the  notation  haB  been  chosen  so  that 
Pi  is  the  source,  Pn  the  sink.  For  certain  values  of  J  -  l...,n 
we  shall  define  labels  v^ , recursively  as  follows. 

Let  Vj  -  oo  ,  A*i  -  0.  For  those  J  ouch  that  a^j  >0,  define 
VJ  *  alJ,/MJ  "  1  *  8®nerai»  from  those  i  which  have  received 

labels  v^,/^  but  which  have  not  previously  been  examined,  select 
an  1  and  scan  for  all  J  such  that  a^j  >  0  and  v y  M  have  not 
been  defined.  For  these  J,  define 

(6)  Vj  -  min  (v1#  a^),/^  *  *• 


Continue  this  process  until  vn, M  have  been  defined,  or  until 

no  further  definitions  may  be  made  and  v  have  not  been 

n  n 

defined.  In  the  latter  case  the  computation  ends.  In  the  former 
case,  proceed  to  obtain  a  new  a^  matrix  as  follows. 

In  general. 


Repl&oe  a.  _  by  a..  — v  and  by  anii 

Ann  Mnn  n  % 


n 


i 
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replace  a^  j  J  ~  vn  aJ/u  aj^  +  vn '  where  ea°h  J 

J  J  J  J 

Is  the  My  of  the  preceding  J'  in  the  backward  replacement. 

This  replacement  continuee  until  /<  -  1  has  been  completed. 

The  labels  Vj,  are  then  recomputed  on  the  basis  of  the  new 

A  matrix  and  the  process  is  repeated. 

Notioe  that  vn  ia  a  positive  integer,  and  hence  the  process 

terminates.  Upon  termination,  the  maximal  flow  X  la  given  by 

defining 

(7)  1  -  max  (ClJ  -  a1J,  0), 


as  we  now  prove 


Lemma  1.  X  Is  a  flow. 

Proof .  Clearly  0  £  x^j  °1J#  Bljnce  a^j  £  0.  It  remains 
to  show  that  X  satisfies  £  (x^  -  x  )  -  O  for  i  -  2,...,n-l. 

J 

Now  the  prooeas  ensures  that  c.  .  +  c..  -  a +  aM.  It, follows 

***  J  J-*-  i-J  Ji. 

from  this  and  the  definition  of  X  that 


(xij  •  X  (0u  au )  • 


It  therefore  suffices  to 


prove  that  ^  is  Invariant,  for 


1  -  2 ,  .  .  . , n— 1 ,  under  the  computation,  as  certainly  ^  (c^  -  a^j ) 
for  the  starting  point  A  -  C.  But  If  1  (/  l,n)  is  the  of 

some  i,  then  there  is  a  k  -My  Thus,  for  this  i,  the  new  a'^ 

are  either  equal  to  the  old  a^j  or  are  given  by 

f alj  -  vn  for  •>  ■ 

alJ  -  <  alj  +  vn  for  J  -  k- 

a*  .  otherwise ; 


* 
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hence 


Lemma  2.  X  It  a  maximal  flow. 

Proof .  At  the  point  where  termination  occurs,  we  have 
defined  a  set  S  of  nodes,  consisting  of  thoso  nodes  for 
which  v^,a^  have  bean  defined;  and  further,  Pi  ^  S,  Pn  <£  S. 
Consider  the  set  p  of  directed  arcs  such  that  3, 

Pj  4-  S. *  Clearly  -  0  for  such  pairs  i,J,  as  otherwise 

we  would  have  defined  v^,  /<  . 

We  will  show  that  p  la  a  cut  whose  value  la  equal  to 
Y  x.  . ,  thus  proving  that  X  is  a  maximal  flow  and  P  a  minimal 


cut . 


That  P  is  a  out  Is  clear.  Por  If  there  were  a  chain 


PiP.  ...P.  P_  with  the  area  P4P,  £  P  ,...,  P,  <$■  P  , 

A.  [  **  ^  l  •* 

then  we  could  deduce  successively  (Blnce  S)  that  P^CS,...,  Pn£S, 

a  contradiction.  To  see  that  P  has  value  equal  to  the  flow 

value ,  notice  that 


l  (cu  -  aijJ  -  1 


(1  <  1  <  n), 

(1-1)  . 


Now  sum  these  equations  over  those  1  for  which  pps.  On  the 
left  side,  If  P1  and  are  both  in  S,  then  c^  -  and 
c  _  a  are  both  In  the  summation  and  are  negatives  of  each 

J  A  J  X 

other.  All  that  remains  are  terms  of  the  form  c^  -  a1 . ,  where 
ppfS,  Pj  <t  S.  Por  these,  as  we  pointed  out  above,  -  0. 


in  [5]. 


TYie  set  pis  actually  the  set  of  "left  arcs"  defined 
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Thus  the  sum  on  tha  loft  reduces  precisely  to  the  sum  of 
capacities  of  members  of  P  ,  as  was  to  be  Bhown. 

We  append  a  simple  example,  using  the  matrix  format  Instead 
of  a  picture  of  the  network,  to  illustrate  the  computation  in 
this  form. 


Example . 


A»  -  C 


As 


As 


A« 


A» 


4  2 

0  4 

0  4 

0  2  2 
0  0  0 


i  *- 
0 

0 


lO 


2 

4 

2 

2 


2 

4 

2 

0 


2  1 
2  4 

1  4 

0  2  2 

0  2  i 

2  1 

2  4 

1  4 

2  2  2 

u0  2  1 


4  4 

1  4 

,2  4  2 
L°  2  1 


2  01 
2  2 
3  1 

4 

0  J 

2  0 
2  0 
3  1 

4 


2  O' 
2  0 

3  0  | 
4 

0 

0  o'1 
2  0 
3  0 
2 

2 

0  O' 
0  0 


V  r* 

ao  0 

4  1 

2  1 

2  1 

2  2 

oo  0 

2  1 

2  1 

2  1 

1  3 

ao  0 

2  1 

1  1 

2  1 

2  4 

oo  0 

2  1 

1  1 

2  2 

2  4 

ao  0 

1  3 

1  1 

i  3 
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Thus  the  computation  terminates  with  Ao,  giving  a  total 
flow  of  7*  IVie  minimal  cut  comprises  the  area  15#  25#  35#  45, 
with  value  7. 

4  .  TUB  HITCHCOCK  PROBLEM 

The  Hitohcock  transportation  problem  la  perhapa  one  of 
the  "most  solved"  linear  programming  problems  in  existence . 

We  shall  propose  yet  another  computation  for  the  problem  which 
will  amount  to  solving  a  sequence  of  flow  problems  of  a  partic¬ 
ularly  simple  kind.  The  basic  idea,  which  stems  from  a  proof 
given  by  Egorvary  [4]  for  a  theorem  of  Konlg  [o ,  p.  232]  on 
linear  graphs,  has  been  used  by  Kuhn  [7]  to  develop  a  very 
efficient  combinatorial  algorithm  Tor  the  optimal  assignment 
problem,  a  speoial  case  of  the  Hitchcock  problem.  Our  method 
differs  only  in  details  from  the  Kuhn  algorithm  in  this  case. 

We  take  the  Hitchcock  problem  in  the  following  form:  Given 

a  matrix  D  -  (d.  .)  of  nonnegative  Integers,  and  two  sots  of 

'  m  n 

nonnegative  integers  (b1#...,b  ),  with  £  ai  "  I  b 

m  n  1-1  x  J-l 

it  is  desired  to  find  a  matrix  X  -  (x^)  satisfying  the  constraints 

X1J  *  °' 

(8)  ^  *iJ  "  &1  * 

l  XiJ  *  bJ 


which  minimizes  the  linear  form 
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A  physical  Interpretation  of  the  problem  Is  that  there  are  m 

t*  h 

originating  points  for  a  commodity,  the  1  point  having  a^ 
units,  and  n  destinations,  the  one  requiring  bj  units.  If 
djj  Is  the  cost,  per  unit  of  commodity,  of  shipping  from  origin 
1  to  destination  J,  find  a  shipping  program  of  minimal  cost. 

Ths  dual  of  the  Hitchcock  problem  la:  Find  cC^,  satis¬ 
fying  the  constraints 

(10)  ^i  -*■  <  d^j  (i  -  l,...,m;  J  -  l,...,n) 

which  maximize 


(11)  I  ai  ♦  I  bj  (?y 

1  J 

The  proof  that  the  algorithm  to  be  described  yields  a  solution 
to  the  Hitchcock  problem  will  be  based  on  the  fact  that  the  dual 
form  (11)  increases  by  at  least  one  unit  with  the  solution  of  each 
successive  flow  problem. 

Eaoh  of  the  flow  problems  will  be  of  the  following  form. 

Find  X  ■*  (x^j)  satisfying 


(12) 


*  0 


I  <; 

J 


^  *> 


j  ' 


-  0  for  a  given  set  _fl-  of  pairs  i,J, 


which  maximizes 


1J* 


(13) 
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This  is  aotually  a  special  Hitchcock  problem.  To  see  that  it 
is  also  a  flow  problem,  set  up  the  directed  network  of  n  +  n  +  2 
nodes  shown  in  ?ig.  6, 


Fig.  6 


Source 


k 


where  the  capaolty  on  the  directed  aro  P^j  is  zero  if  , 

large  otherwise,  the  capacities  of  source  and  sink  aros  are  the 
and  bj,  as  shown,  and  interpret  x^j  as  the  flow  from  to  Qj . 

To  solve  such  a  problem  we  may  of  course  use  the  computation 
of  Section  5  involving,  in  this  case,  an  b  +  n  +  2  by  m  +  n  +  2 
matrix.  It  is  possible  (and  computationally  convenient)  to  de¬ 
scribe  the  process  in  terms  of  an  m  by  n  array.  The  verification 
that  the  two  descriptions  are  the  same  for  this  particular  clasa 
of  problems  will  be  left  to  the  reader. 

Let  X  be  a  solution  of  the  constraints  (12).  Corresponding 
to  certain  of  the  rows  1  -  l,...,m  of  X  we  will  define  integers 
v1,/^lj  similarly,  for  certain  of  the  columns  J  •  l,...,n  we  will 


•74  3 
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deflne  Integers  wy  *y  These  definitions  will  be  made  re¬ 
cursively;  first  a  set  of  vlf  ^  will  be  defined,  from  these  a 
set  of  w.,  A 1  will  be  defined,  and  so  forth,  alternating  between 

the  rows  and  columns . 

For  those  i  such  that  £  <  al»  define 

(14)  vi  "  ai  "  I  xiJ*  A1  "  °  ‘ 

J 

Mext  select  an  1  which  has  been  labeled  and  scan  for  all  (unlabeled) 
j  such  that  l,J<fcn;  for  these  J,  define 

(13)  -j  -  V  ^  -  i. 

Repeat  until  the  prevloualy  labeled  1'.  are  exhausted.  We  then 
select  a  labeled  J  and  scan  for  unlabeled  i  such  that  Xjj  >  0; 
for  these  1,  define 

(l6)  -  min  (x^,  "jK^i  " 

Repeat  until  the  previously  labeled  J • •  are  exhausted.  Again 
select  one  of  the  i's  just  labeled,  look  for  unlabeled  J  such 
that  ifj*Jl  ,  and  define  Wj ,  Aj  by  (15).  Continue  In  this  fashion, 
using  (15)  and  (l6)  alternately  until  either  Aj  have  been 
defined  for  some  j  with  7  x^  <  bJ(  or  until  no  further  defini¬ 
tions  may  be  made.  In  the  latter  case  X  is  maximal;  in  the 
former  we  can  get  an  improvement  aa  follows.  Let 

(17)  v  -  min  (Wj,  bj  -  ^  *ij)  * 

Alternately  add  and  subtract  v  from  the  sequence 


f 


’-7*3 
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(18) 


V  XajJi#  XliJi'  "*,Xlk-iJk,JClkJic  ' 


wiis ro  J  |  ® ^  i  1 1  *  .  p  j •  ■  ^ ®*  /S  .  |  *  *  *  i  J i.  *  j 

A  j  J »  1  *  J  ■  *  xk-l  ' 


lk  -  Aj  ,  and  -  O . 
Example  . 
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Celia  containing  a  circle  comprise  -TL  (the  complement  of 
jTL  ).  The  entries  within  the  clroles  (zeros  elsewhere)  con¬ 
stitute  an  X  satisfying  the  constraints  (12).  The  defining  pro¬ 
cess  terminates  with  w*  -  1,  >*»  -  4,  since  column  5  1b  a  column 
in  which  the  sum  of  the  entries  is  less  than  b^;  l.e.,  3  +  5  <9* 


P-7^3 


1 2-2  9^5 


The  sequence  along  which  an  Improvement  of  min  (wB,ba  —  V  x^  )  "  1 
can  be  made  la  x4a  -  x4l  -  1,  xTi  -  0,  xTs  -  1,  x.B  -  0,  ao 
la  easily  read  off  ualng  the  ^'s  and  A^a  alternately. 

We  are  now  In  a  position  to  describe  a  general  routine  for 
the  Hitchcock  problem.  Aa  a  starting  point,  form  the  difference 
matrix 


(19) 


du  -  ^  • 


where  -  min  d^j,  "  min  (d^  -  Qt .  Ttius  d^  —  ^  0; 

l*e.»  (X^,  0^  8*tlafy  the  dual  constraints  (10).  The  next  step 
Is  to  solve  the  flow  problem  with 


(20) 


J\  -  [iJ|d1J  -  a1  -  pj  >  0 


If  the  maximizing  flow  X  satisfies 


i.J 


tjj  •  ^  ai*  then  X  la  a 


minimizing  solution  to  the  original  Hitchcock  problem.®  If,  on  the 

8  Hila  la  easily  deduced  as  follows.  Obaerve  first  of  all 
that  for  any  c^,  0  ,  the  two  Hitchcock  problems  with  cost  matrices 
dlJ  0X1(1  dij  ~  ^i  ~  are  equivalent;  for  ^  x ^  -  a^ 

\  hj  -  bj  ^  th»*  £  <au  -«!  -  *u  -  ,1/u-u 

—  al  ~  ^  b  |  ^ j »  ^d  the  last  two  sums  on  the  right  are 

Independent  of  x^ . .  Thus  If  we  have  oh,  f  with  d.  —  d,  —  >  0, 

and  are  able  to  find  an  X  satisfying  (8)  and  x..  -  0  for  1J£/L  , 

^  J 

then  clearly  X  minimizes  £  (d^j  —  hence  minimizes 

Zi  -  1,J 


u  xi  y 


P-7.^3 
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other  hand,  ^  X±J  <  £  a±,  lot  I  bo  the  Index  set  of  the  labeled 

rows  of  X,  J  the  index  set  of  the  labeled  columns,  and  define  new 
dual  varlablea  by 


+  k 


(21) 


"k 


(1  *  I), 

U  i  I)  , 

(J  c  J), 

(J  4  «J )  # 


where  k  -  rain  (d^  -  ^ ) .  Notice  that  k  >  0,  Blnoe  pairs 

1#J  with  ifcl,  are  contained  in  _Tl 
LSHSi  J.  !«!«;♦  £  bj*'  >f 

Proof,  Hie  fact  that  X  la  a  maximal  flow  implies  that 
^  X1J  *  for  Jejr  •  Also,  the  labeling  process  ensures  that 
X1J  "  ®  ^or  *4 1#  and,  as  we  have  mentioned  earlier,  x^  -  0 

for  1£I,  J4J.  Since  all  i  with  7  x ^  <  a±  are  in  I,  it  follows 
that  J 

lii  >  xh  :h  ■  jL  bJ  • 

and  hence  the  new  dual  form  haa  been  increased  by  the  amount 


1 A  ■>  -  1  >  °  • 


I 
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Another  way  to  see  this  Is  to  note  that  the  minimal  cut  in 

the  associated  network  (Fig.  6)  has  value  ^^ai  +  ^  ^  ai* 

and  hence  T  a.  >  T  b ,  . 

i Cl  1  jtj  J 

Now  fora  the  new  difference  matrix  d^j  —  cC^  —  /?j  >  0  by 
subtracting  k  from  the  I— rows  of  the  previous  difference  matrix, 
and  adding  k  to  the  J— columns .  We  may  now  take  the  maximal  flow 
X  of  the  previous  flow  problem  as  a  starting  point  in  the  new 
flow  problem  and  proceed  as  before  . 
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